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Abstract 

We show that A-symmetries can be algorithmically obtained by using the Jacobi 
last muItipUer. Several examples are provided. 
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1 Introduction 

Lie group analysis is the most powerful general tool to find the solution of Ordinary 
Differential Equations (ODEs). However an ODE of nth order does not always admit 
Lie point symmetries. Moreover this Lie group analysis is useless when applied to n 
equations of first order because they admit an infinite number of symmetries, and there 
is no systematic way to find even a one-dimensional Lie symmetry algebra. One may 
try to derive an admitted n-dimensional solvable Lie symmetry algebra by making an 
ansatz on the form of its generators but when successful (rarely) it is just a lucky guess. 
However, in [26] it has been shown that any system of n equations of first order can be 
transformed into an equivalent system where at least one of the equations is of second 
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order. Then the admitted Lie symmetry algebra is no longer infinite-dimensional, and 
nontrivial symmetries of the original system may be retrieved. This idea has been 
successfully apphed in several instances ([26], [36], [31], [27], [16], [28], [5], [8]). 
In [29] another method was devised. It uses the Jacobi last multiplier ([11], [12], [13], 
[14], [17], [18], [1], [37]) to transform a system of n first-order equations into an equivalent 
system of n equations where one of the equations is of second order, namely the order of 
the system is raised by one. In [29], among other examples, the method was successfully 
applied to the second-order equation [15] [Ch. 6, 542j9^: 



where p 7^ is a real constant and / 7^ is an arbitrary function of the independent 

variable t. This equation docs not possess Lie point symmetries for general f{t) and 
yet is trivially intcgrable [7] . In [29] the introduction of the Jacobi last multiplier led to 
an equivalent system of two equations, one of first order and the other of second order, 
which admits enough Lie symmetries in order to integrate it by quadrature. 
In [20] and [21] Muriel and Romero introduced the so-called A-symmetries which were 
later included into the telescopic symmetries by Pucci and Saccomandi [35]. Again as 
in the case of Lie point symmetries of first-order equations the real problem is to find 
solutions of the determining equations. This problem has been clearly stated by Pucci 
and Saccomandi on page 6154 of their article. 

In [4] and [6] Cicogna, Gaeta and Morando provided a geometrical characterization of 
A-symmetries for both ordinary and partial differential equations. However, in either 
paper they did not address the problem of solving the determining equation, they just 
obtain a A-symmetry by guesswork. 

In [2] Catalano Ferraioli showed that A-symmetries correspond to a special type of non- 
local symmetries and in a final remark observed that also telescopic symmetries can be 
recovered by nonlocal symmetries. Yet there is much guesswork involved in order to find 
nonlocal symmetries, and on page 5485, Remark 2, Catalano Ferraioli wrote: "However, 
we note that the problem of determining the general solution of (11)^, should be at least 
as difficult as solving the given ODE. Therefore in practice it could not be so easy to 
determine such correspondence (see example 3)" . 

In [30] the method described in [29] was applied to example 3 in [2] , i.e 



that is a second-order Painleve-type equation [34] , which does not admit any Lie point 
symmetries and is a particular case^ of Painleve XIV equation [9], i.e.: 



y"^— + nt)y^^' + pf{t)yV, 



(1) 




(2) 




(3) 



^ It is a linear first-order partial differential equation. 
^It corresponds to assume Q{t) = 1 and S{t) = t. 
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which does not possess any Lie point symmetry for arbitrary Q{t), S{t) although it has 
a Riccati-type first integral, i.e.: 

y' - Q{t)y^ + S{t) 

= (4) 

with Oi an arbitrary constant. In [30] the introduction of the Jacobi last multiplier into 
equation (2) led to a system of three equations, two of first order and one of second order. 
This system admits a three-dimensional solvable Lie symmetry algebra and therefore can 
be reduced to a Riccati equation that can be integrated in terms of Airy functions. Thus 
a new first integral different from (4) was found, i.e. 

(^2 ^y,^ t)kiryM (0 + 2yAiryAi (1, Q + 2{t - y')y^ + {t + y'f 

(y2 + y' + i)AiryBi(0 + 2yAiryBi(l,0' ^ V 

(5) 

with a2 an arbitrary constant. Then, combining the two first integrals (4)- with Q{t) = \ 
and S{t) = t - and (5), one gets the general solution of equation (2) in implicit form, 

{2y + gQAiryAi (Q + 2AiryAi (1, Q ^ A = f + ^ ffi^ 

{2y + ai)AiryBi (0 + 2AiryBi (1, 4 ' ^ ' 

This suggests that there might be a direct link between the Jacobi last multiplier and 
A-symmetries. 

In the present paper we show that the introduction of the Jacobi last multiplier 
allows to find A-symmetries algorithmically. 

In the next Section, we briefly recall the connection between nonlocal symmetries 
and A-symmetries as shown by Catalano Ferraioli in [2], and some essential properties 
of the Jacobi last multiplier. Then we show how to use the Jacobi last multiplier to find 
A-symmetries. Recently in [23] Muriel and Romero have classified all the A-symmetries 
of any second-order ODE through an equivalence relationship and proved that two A- 
symmetries lead to functionally independent first integral if and only if they are in 
different equivalence classes. This important result allows us to discriminate the more 
than one A-symmetry that we found in some of the examples. In Section 3 we apply 
our method based on the Jacobi last multiplier and either recover known or find new 
A-symmetries of several differential equations of second order: equation (1) as given in 
[7], equations Painleve V, XIV (3), XV, and XVI as given in [9]^, and examples 4 and 5 
as given in [3]. Section 4 is devoted to some conclusions. 

2 A-symmetries and Jacobi Last Multiplier 

Let us consider an n*'*-order ODE: 



^In particular we recover the known first integrals [9]. 
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where by an apex we mean the order of differentiation. In [2] a nonlocal interpretation of 
A-symmetries was given. There Catalano Ferraioli has shown that seeking A-symmetries 
of eq. (7) is equivalent to add to equation (7) the equation 



A 



(8) 



for the new field u{t) with A = X{t,y,y\y", . . . ^)). The symmetries of (7, 8) are 
obtained by considering the infinitesimal operator: 



Y = f{t, y, u)dt + fi{t, y, u)dy + i{t, y, y', y", y^''-^\u)d^ 
with the constraint 

Under the constraint (10) we get Y — e'^[X -\- ^dj\ where 

X = T{t,y)dt + 'n{t,y)dy 
and ^ = ^{t, y, y', y", . . . , y^'^~^^). The prolongation of X is 



(9) 
(10) 

(11) 



prX = T{t,y)dt + r]{t,y)dy + r]^'\t,y,y',y" 
+V^'\t,y,y',y",...,y^^-'^)dy. + .. 



,("-!) 



)9y' 



with 



(A + A)r/(")-i/'(A + A)r 



(12) 



(13) 



where Dt = dt + J2k=oy^''^^^'^v'^''^' y^^'' ^ ^' = rj. A A-symmetry of equation (7) 

is any solution of the determining equation 



prX (y W - fit, y, y', y", y^^''^) = 0, 



(14) 



depending on the three unknowns A, r and rj and therefore is highly undetermined. In 
[23] proved that any A-symmetry can be put in an evolution form as a translation in y, 
i.e. with T — and rj — 1. Hence equation (14) becomes determined since it has only 
one unknown, but it is still difficult to solve it since it is a nonlinear partial differential 
equation in the unknown A. 

If equation (7) is transformed into an equivalent system of first-order equations, i.e. 

= Wj(t,wi,...,w„), (15) 
then its Jacobi last multiplier M is obtained by solving the following differential equation 



dlog(Af) 
dt 



i=l 



dwi 



(16) 



I.e. 



M = exp 




(17) 



In [29] all the properties of the Jacobi last multiplier are listed. There three strategies 
were proposed with the purpose of finding Lie symmetries of any system (15): 

1. Eliminate - if possible - each of the variables Wi in order to obtain an equivalent n- 
order system which contains a single equation of second order and n—2 equations of 
first order. The admitted Lie symmetry algebra is no longer infinite-dimensional 
and the Lie group analysis can be usefully applied ([26], [36], [5]). Prom this 
strategy we also get first integrals ([19], [27], [16]). 

2. Decrease the order of system (15) by one choosing one of the variables Wi as the 
new independent variable. Then apply cither Strategy 1 or 3 in this list. For 
example, if Wi = y is the new independent variable, then system (15) becomes 



This method has been applied in several examples (e.g. [31], [16], [28]) since its 
first instance, the Kepler problem [26] . 

3. Increase the order using the transformation suggested by the Jacobi last multiplier, 
i.e. introducing a new dependent variable R such that 



and eliminating - if possible - each Wi which appears in (19). Then system (15) 
reduces to a system given by a single second-order ODE and n — 1 first-order 
equations. If a new independent variable was chosen (Strategy 2), say wi = y, 
then R will satisfy the following equation: 



dwk Wk 



^k{y,W2, ■■■,Wn) 



{k = 2,n) 



(18) 



dy W, 




(19) 




(20) 



and then the strategy goes as above. 
In [29] Strategy 3 was apphed to find Lie symmetries of several systems. 



The similarity between Strategy 3 and the nonlocal approach to A-symmetries as given 
by Catalano Ferraioli suggests to search for A-symmetries such that 



This implies, when feasible, that oo = log(l/M). In fact this connection cannot be made 
if the divergence of the system (15), namely Div = Yl^=i zero, since then any 

Jacobi last multiplier is a first integral of (15) and therefore such is ou. 
However there are many ways in which a system of n first-order equations can be written 
as a single equation of nth order and viceversa. The Jacobi last multiplier is then 
different and if one way yields Div = 0, another way may yield Div ^ 0. In particular 
the following system of two first order equations 



^ = W^2(i,wi) 



(22) 



has Div — and therefore M[^^ ,„2] = 1 is one of its Jacobi last multipliers. If we derive 
w\ from the second equation, i.e. w\ — Wiit^ Wg), then an equivalent second-order ODE 
is obtained, i.e. 



which has Div ^ since^ 

M - M 9{Wi,W2) 
^ ' ^ 'o{W2,W2) 



dt 



d 



dw 



7^2(t,^) 



(23) 



^jW2{t,w'2) 



dw. 







d 



dw. 



jW2{t,w'^). (24) 



i.e. its Jacobi last multiplier cannot be a constant. An illustrative example of such an 
instance is the following system studied in [33] 



6exp(r2) + a 



r'^ — B exp(ri) + A 



(25) 



which has obviously Div — 0. Following [33] we can transform system (25) into an 
equivalent second-order ordinary differential equation by eliminating ri. In fact from 
the second equation in (25) one gets 



ri = log 



B 



(26) 



^This is one of the known properties of the Jacobi last multipher [1], [29]: given a non-singular 
transformation of variables 

T : {wi,W2,. ■ ■ ,Wn) — > (ri,r2, . . . ,r„), 
then the last multiplier Mj^j of the new system r'^ = Ri{t, ri, . . . , r„) is given by: 



Mr, 



Mr, 



d{ri,r2, . . . ,r„) 
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and the equivalent second-order equation in r2 is the following 

r^' = -(6exp(r2) + a)(A-r^). (27) 
which has Div = 6exp(r2) + a 7^ 0. 

As a final remark we recall that recently [23] Muriel and Romero have proved the equiv- 
alence between two A-symmetries of a second-order ordinary differential equation 

y"^<f>{t,y,y'). (28) 

Assuming that Ai and A2 yield two A-symmetries Xi = Ti{t,y)dt + rii{t,y)dy and X2 = 
T2{t,y)dt + ri2{t,y)dy of equation (28), respectively, they are equivalent if and only if 

QiiA + X2){Q2) - Q2{A + Ai)(gi) = 0, (29) 

where Qi = fji — y'ri, {i = 1,2) and A = dt + y'dy + (f){t,y,y')dyi is the vector field 
associated with (28). This equivalence will be used in the next Section 3 when we 
obtain more than one A-symmetry. 



3 Examples 

In the following examples we denote by A^ the A-symmetries presented by Muriel, 
Romero, Catalano Ferraioli and Morando in the references [20], [21], [3] and by Aj 
the A-symmetry that we find by using formula (21) once we rewrite (7) as (15). 

3.1 Equation (1) 

In [20] a A-symmetry of equation (1) was determined, i.e. 

X^>^''^^dy With Xk^py''f{t)+y'/y. (30) 
The divergence of equation (1) yields 

Xj^py^f{t) + 2^. (31) 

If we put A J into (12) then the solution of the determining equations (14) yields two 
A-symmetries, i.e. 

Xi'^ = -dy, ) = \dt + y^-'mdy. (32) 

y y 

The first prolongation of X^\ i.e. 

prX?^ = X}^) + (py'-'fit) + ^) dy, (33) 



yields the first-order invariants 



= + ^1 = ^ (34) 
y 

that replaced into equation (1) generate the first-order equation 

as obtained in [20]. The first prolongation of ^2'^^ i.e. 

prX(") = + {y'^-'fit)p + y^-'f\t) + y'-'y'fit)) dy, (36) 
yields the first-order invariants 

where F{t) = t2—p J f{t) dt, that replaced into equation (1) generate the same invariant 
as Xi'l 

This result should not be a surprise since the two A-symmetries (32) and the A-symmetry 
Afc found in [20] are all equivalent as defined by Muriel and Romero in [23]. In fact 
substituting x[^^ and Xg^^ - that have the same Aj (31) - into (29) yields: 

QM + A,;)(Q2) - Q2{A + \j){Qi) = Qi{A){Q,) - Q2{A)iQ^) 

1 " 



y 



y'-'f{t) + y' ((p - i)/-V(^) + 2^) - ^ (7 + f'^^^y'^' + pf{t)yy) 

+ {y'-'m - ^) ^ - (38) 

since Qi = 1/y and Q2 = y^^^f{t) — y' /y"^. Also substituting x'^'' and X^'^'^^ and their 
corresponding Aj and A^ into (29) yields: 

Qi{A + \u)m-Qk{A + \j){Q,) = -- (py^fit) + y^]+^(pyPf{t) + 2y^)-^^= 

y \ y J y\ y J y^ 

(39) 

since Qk — 1- 

3.2 Painleve-Ince V 

The divergence of Painleve-Ince equation V 

y"^-2yy' + q{t)y' + q\t)y (40) 
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yields 

\j = -2y + q{t). (41) 
If we put \j into the A-prolongation then (14) yields two A-symmetries, i.e. 

= dy, = 9, + {yq{t) - dy. (42) 

These two A-symmetries (42) are equivalent [23]. In fact 

QM){Q2) - Q2{A){Q,) = yq'{t) + y' {q{t) - 2y) - {-2yy' + q{t)y' + q'{t)y) = (43) 
since Qi = l and Q2 = yq{t) - 2/^ - y'. 

Therefore we consider only one A-symmetry, x['^\ Its first prolongation, i.e. 

prXf ) = + {q{t) - 2y) dy, (44) 
yields the first-order invariants 

yi = -y(i{t) + + y\ h = t (45) 

that replaced into equation (40) generate the first-order equation 

y[^O^y^^ai^ -yq{t)+y^ + y' ^ai, (46) 

and therefore the known first integral of (40) is derived [9] . 
As far as we know x[^^ in (42) is a new A-symmetry of (40) . 

3.3 Painleve-Ince XIV 

In [22] a A-symmetry of equation Painleve-Ince XIV equation (3) was determined by 
assuming that A was linear with respect to y', and thus the following A-symmetry was 
found 

= with Afc = yg(t) + ^ + ^. (47) 

y y 

Instead the divergence of equation (3) yields 

Xj^^ + Q{t)y + Dt{\og{y')). (48) 

y 

If we put A J into the A-prolongation and solve the determining equation (14) we get 
two A-symmetries, i.e. 

Xi'^ = '-dy, ) =. + + Qit)^ dy. (49) 
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One can prove that these two A-symmetries and that found in [22] are equivalent. 
The first prolongation of x['^\ i.e. 

yields the first-order invariants 

Sit) , y' 
yi = ^-Q{t)y + -, h = t 
y y 

that replaced into equation (3) generate the first-order equation 

= =^ yi = oi =^ Q{t)y + — ^ai, 

y y 

and thus the known first integral (4) of (3) is derived. 

3.4 Painleve-Ince XV 

Painleve-Ince XV equation 

// y'^ y' / N 2 d fr'it) 



y y dt\ r{t) 

is known to possess a first integral [9], i.e.: 
1 fr'{t) 



^-y + y' + lj - 2 (^r(t)y + j r(t)dtj 



ai. 



1/2 \ r{t) 

The divergence of Painleve-Ince XV equation (53) yields 

A^ = - + A(log {y')). 
y 

If we put Aj into (12) then (14) yields one A-symmetry, i.e. 

y' ' rW 

The first prolongation of X^^\ i.e. 

prX(^) = X(^) + -^'(^)^(^)^' + ^'i^fy^ - r'{t)r{t)y{y' + 1) - r{t)\y' + 1) ^ ^ 

yields the first-order invariants 



^ r'(t) y' + l ~ , , /• , , , 

y = 4-( + ^ , t = r{t)y + / r{t)dt 

r[t) y J 



10 



that replaced into equation (53) generate the first-order equation 

I = i ^ - 2f = ^ + t±iy - 2 (r(t)y + / r^*) = a„ (59) 

and thus the known first integral (54) of (53) is derived. 

As far as we know X^^^ is a novel A-symmetry of equation (60). 



3.5 Painleve-Ince XVI 

Painleve-Ince XVI equation 

y" = - - q'it)^ + y'- q{t)y' + q{t) (60) 
y y 

is known to possess a first integral [9], i.e.: 

= (61) 
The divergence of Painleve-Ince XVI equation (60) yields 



A. = -^ + A(log(y^)). (62) 



If we put Aj into (12) then (14) yields one A-symmetry, i.e. 
The first prolongation of X^^\ (12), i.e. 



yields the first-order invariants 



that replaced into equation (60) generate the first-order equation 



(63) 



f = ?/-g(t) (65) 

y 



- -^y -t =ai^ -{y-q{t)) =ai, (66) 



dt y \ y J 

and thus the known first integral (61) of (60) is derived. As far as we know X'^^^ is a 
novel A-symmetry of equation (60). 
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3.6 Example 4 in [3] 

In [3] a A-symmetry of equation 

y" = ity' - ty" + y") exp(-l/|/) + 2y + (67) 

was determined with = texp(— 1/y) — 1/t. This equation has been completely solved 
by considering its solvable structures [3]. 

The divergence of equation (67)- that has no point symmetries - yields 

\j = t exp(-l/y) + A (log (y') + t) . (68) 

If we put A J into the A-prolongation and solve the determining equations (14) we get 
two A-symmetries, i.e. 

'(A) _ 1 a 1 / 1 . t 



' y^eMt) ' ' y^ewm W exp(l/y) V 

One can prove that these two A-symmetries and that found in [3] are equivalent. 
The first prolongation of x'^\ i.e. 

' ^ \y^exp{t) y^exp{t + l/y)J " ^ ' 

yields the first-order invariants 

7/i = ^-texp(-l/|/), h = t (71) 

y 

that replaced into equation (67) generate the first-order equation 

= yi =^ yi = aiexp(i) - iexp(-l/y)^ exp(-i) = oi, (72) 

and therefore a first integral of (67) is derived. 

3.7 Example 5 in [3] 

In [3] a A-symmetry of equation 

y" = 2^ + [teMt/y) --^y'- + 1/) exp(t/y) + ty' + 2^ (73) 

was determined with A^ = texp{t/y). The divergence of equation (73)- that has no 
point symmetries - yields 

A^ = texp(t/y) + A (log (^^))) . (74) 
12 



If we put Xj into (12) then (14) yields two A-symmetries, i.e. 

One can prove that these two A-symmetries and that found in [3] are equivalent. 
The first prolongation of x['^\ i.e. 

prXf") = ) + ^ {t^y exp{t/y) -y + 2ty') dy, (76) 

y 

yields the first-order invariants 

ty'^exp{t/y) -y + ty' 
ty 

that replaced into equation (73) generate the first-order equation 



yi = , = t (77) 



, -3yi + ^ 5ai + ^ ,2 ^2/' exp(Vy) - y + ^2/' 
= ^ ^y^^^^^^ f 5 ^^^^ 

and therefore a first integral of (73) is derived. 

4 Conclusions 

In this paper we have shown that the Jacobi last multiplier provide and algorithmic 
simple way to construct A-symmetries. Once a A-symmetry is obtained then it is a 
simple task to derive a first integral as we have shown in the many examples presented 
in Section 3. 



We remark that Strategy 3 as described in [29] is not equivalent to the reduction by 
using A-symmetries. In fact different output are obtained. For example, equation (1) 
was completely solved by quadrature using Strategy 3, while the reduction by using 
A-symmetries yields just a first integral. Yet those two methods may complement each 
other as in the example of equation (2). 



Also one should be aware of the fact that a A-symmetry could be equivalent to a Lie 
point symmetry. Lie point symmetries may be obviously considered A-symmetries with 
A = 0. For example, let us consider equation (38) in [22], i.e. 

2yy"-6y" + y'' + y' = 0, (79) 

which admits a trivial Lie point symmetry, i.e. T — dt- The divergence of this equation 
yields 

Xj = 6^ = A(log(/)). (80) 
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If we put Xj into (12) then (14) yields one A-symmetry, i.e. 

X^'^ = \ dt. (81) 

yb 

It is easy to prove that this A-symmetry is equivalent to the Lie point symmetry F. 
We now search for A-symmetrics of equation (27) which, ii A = a, admits two Lie point 
symmetries, i.e. Fi = dt, F2 = cxp{—at){dt + aS^j)- The divergence of equation (27) 
yields 

Aj = 6exp(r2) + a. (82) 
If we put Aj into (12) then (14) yields one A-symmetry, i.e. 

=exp(-at)(9^2. (83) 

It is easy to prove that this A-symmetry is equivalent to the Lie point symmetry F2. 

A more detailed analysis of the Jacobi last multiplier approach to A-symmetries is needed. 
In particular the analysis of higher order ODEs can provide new ideas and confirm the 
importance of this method for the integration of ODEs. Moreover its extension to PDEs 
can provide new insights on the meaning of /^-symmetries [4] , [6] . 

Work is in progress in these directions. 
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